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Abstract 

A transfer matrix formalism applicable to certain reparametrization invariant the- 
ories, including quantum gravity, is proposed. In this formulation it is found that 
every stationary state in quantum gravity satisfies a Wheeler-DeWitt equation, but 
each with a different value of the Planck mass; the value Trip lanck turns out to be 
proportional to the eigenvalue of the evolution operator. As a consequence, the fact 
that the Universe is non-stationary implies that it is not in an eigenstate of Newton's 
constant. 
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1 Introduction 



In non-parametrized Hamiltonian mechanics, there is a set of conjugate variables 
{q l ,Pi} and a time parameter t. Classical trajectories are fixed by specifying initial 
{go} and final {(?}} coordinates, and also a time lapse At. In parametrized mechanics 
there is also a set of conjugate variables {q^jP^} and an evolution parameter r. 
The difference is that in parametrized theories, a classical trajectory is fixed by the 
initial {q$} and final {q^} coordinates alone; the values of r that happen to be 
associated with those coordinates are irrelevant. An example is the case of a free 
relativistic particle, where specifying {x^} at the initial and final points determines 
the trajectory. A field-theory example is general relativity, where the initial and 
final three-manifold is sufficient, together with Einstein's equations, to determine the 
four-manifold between them. Since the {q^} of parametrized theories contain time 
information, applying standard quantization prescriptions is, in effect, "quantizing 
time". This poses no problem when the Hamiltonian of the parametrized theory is 
parabolic, as in parametrized non-relativistic mechanics, or parametrized relativistic 
scalar field theory. But for parametrized theories with hyperbolic Hamiltonians, such 
as the free relativistic particle, or quantum gravity, standard quantization procedures 
can lead to serious difficulties in identifying an appropriate evolution parameter, and 
a conserved non-negative norm [|l|. In some cases, e.g. a free relativistic particle in 
flat space, these problems can be easily overcome; in others, such a free relativistic 
particle moving in an arbitrary curved background, they are much more problematic. 
In quantum gravity, these difficulties are known as the "problem of time." 

In this article I will propose a transfer matrix formalism for quantum gravity, 
and certain other parametrized theories of the type described above. Since a transfer 
matrix, by definition, is an evolution operator, this proposal is intended as a possible 
resolution of the time problem in quantum gravity. The method treats all dynamical 
degrees of freedom of parametrized theories on an equal footing: all are operators, 
none in particular is an evolution parameter. This method will first be illustrated for 
a one-dimensional "universe", whose action is that of a free relativistic particle. The 
formalism will then be extended to minisuperspace-type actions, typical of quantum 
cosmology, and finally (with caveats regarding operator-ordering and regularization) 
to full quantum gravity. 
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2 A One-dimensional Universe 



A free relativistic particle is the simplest example of a system with a reparametrization 
invariant Lagrangian; for this reason it is often used as a "warmup" exercise |], ^ for 
higher dimensional reparametrization invariant theories, such as strings and quantum 
gravity. To briefly recall some of the familiar analogies: The action of a relativistic 
particle Q 



5 = - mo /T v ?? (1) 

is invariant under reparametrizations r — ► /(t); the Einstein-Hilbert action is invari- 
ant under general coordinate transformations x^ x ,fl (x). In phase space, one can 
also write for the particle theory 



S = Jdr(p^-NH) 

h = T^-yp^+o ( 2 ) 

and this leads to the Hamiltonian constraint H = 0, which is just the mass-shell 
constraint. The analogous steps, in the ADM decomposition for gravity, lead to the 
superHamiltonian and supermomentum constraints. Upon quantization, the Hamil- 
tonian constraint becomes a constraint on physical states Hip = 0, which is the 
Klein-Gordon equation in the case of a particle, and the Wheeler-DeWitt equation 
for quantum gravity. 

If one slices spacetime into a series of spacelike hypersurfaces £*r with a parameter 
T = T(x) (the time coordinate) labeling each hypersurface, the following norm is 
preserved by the Klein-Gordon equation: 

< V# >t= / dE^pdrf - iPd^iP*) (3) 

This norm is not positive definite in general (although in flat spacetime one can restrict 
to positive frequency states, in which case the norm is positive definite) .f\ A similar 
construction can be made in quantum gravity. It is possible to slice superspace into 
hypersurfaces such that a norm analogous to (|3]) is independent of the hypersurface, 
Signature convention r/ — diag[— 1, 1, 1, 1]. 

5 More generally, it is possible to identify a conserved non-negative norm in spacetimes with a 
timelike Killing vector field. 
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providing \1/ satisfies the Wheeler-DeWitt equation. It is not clear how to make such a 
norm non-negative. Nevertheless, the analogy to the relativistic particle case suggests 
(and the idea goes back to the classic paper of DeWitt ||) that some coordinate in 
superspace can be interpreted as a time evolution parameter in quantum gravity. 
There have been many subsequent variations and extensions of this idea; a recent 
comprehensive review is found in 

On the other hand, there is one aspect of the relativistic particle example which 
seems quite different from the situation in 4D gravity. In the case of the relativistic 
particle, the observer is obviously external to the particle. The information contained 
in the Klein-Gordon wavefunction refers to measurements that can be made by such 
external observers, who are free to measure the particle position on any spacelike 
hypersurface. In contrast, observers in 4D gravity (at least, human observers) are 
internal to the system in question (the Universe), since we live in spacetime rather 
than superspace. Is it possible to find and quantize a one-dimensional model which 
reflects also this aspect of the 4D case? 

Such a model is obtained by simply reinterpreting the action of eq. ([!]). Consider 
an observer living in a one-dimensional universe parametrized by a single coordinate 
r, who is able to make observations on a 4-component "field" x^(t) in that one- 
dimensional spacetime. The reparametrization invariant action of this field is taken 
to be that of eq. ([!]). The Hilbert space consists of wavefunctions 0(x M ) with norm 

< 0|0 >= J rf 4 x/i(x)0*(x /i )0(x /i ) (4) 

where fi{x) is a measure to be determined. Note that since all components x^ are 
observable, the integral is taken over the full four- dimensional "field" space. The 
fact that all components are observable and all are integrated over in the norm 
(||) is the main (and crucial) difference between the "one- dimensional universe", in 
which none of the observables is to be interpreted as an evolution parameter, and the 
relativistic particle. 

In classical relativistic mechanics, the dynamics of a particle moving in curved 
space, or in some external potential, can be described covariantly by a parametrized 
trajectory x^(t). In quantum physics the configuration space is Hilbert space, and 
dynamics can be presumably also be described as a parametrized trajectory ip^x^, r) 
in Hilbert space. Let us postulate a corresponding Schrodinger equation 

ihd T ip(x,T) = Htp(x,r) (5) 
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where H is an operator, hermitian in the measure fi(x) and invariant w.r.t Lorentz 
transformations, such that 

x^{t)=<^x,t)\^(x,t)> (6) 

is a parametrized solution of the classical equations of motion. In this way, the 
parameter r running along a trajectory of Hilbert space can be identified with the 
variable r parametrizing a certain classical trajectory, corresponding to the motion 
of the center of the wavepacket. Of course reparametrizations of r have no physical 
significance, and the Schrodinger equation above can be made to look covariant w.r.t 
reparametrizations of r by intoducing an "einbein" for the trajectory 

ihd T ^(x,r) = e(r)Hip(x,T) (7) 

However, assuming that H is r-independent, the trajectory through Hilbert space 
depends only on an initial state, and not on e(r), which can always be set to e = 1 
by a reparametrization. From here on we set e = 1. 

I will now show how an evolution operator H and measure fi(x), with the required 
properties, may be obtained from a transfer matrix formalism. The starting point 
is the Euclidean action corresponding to eq. ([[]), obtained rotating the field space 
metric (rj^) from Lorentzian to Euclidean signature, and extracting an overall factor 
of i; i.e. 

iSL oren tz[gfj,u = Vnv] * S Euclid[9 iiv = 3[iv} (8) 

Then, generalizing from D = 4 to arbitrary D, the transfer matrix T e is defined by 



i/;(y,T + e) 



T e ip(y,r) = exp[-H t e/h}ij(y,T) 
d D x fj,(x) exp(— S[y, x]/Veh)ij;(x, r) 



(9) 



where S[y,x] is the Euclidean action of the classical trajectory connecting (and ter- 
minating at) points x M and y^ . Since S[x,y] = S[y,x], and S[y,x] is real, it follows 
that H e is hermitian in the measure \x. Note that by the usual trick of integration by 
parts, we have 



< 



Y[d D x fi(x n )j | 

f 1 dS[{ Xi }] 
IK dxt 



dS\{Xi\] 8 , .1 r orr . i—r 

- -Tyi In /M{x h ) j exp [-6 [{Xi}\/ V eh 



eh dx 11 



d 
dxj 



ln/i(x fc ) 



> 



(10) 
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where 

dS[{x, t }] d 



;)../' ~ f)^ ^2 S[ x n+1, x n] 

ux k ux k n 



SS[x(t)) 



(11) 



5x»(t') 

and where 5*[x(r)] is the continuum action of the path x^{r). Formally then, in the 
e — ► limit, we obtain the quantum version of the Euclidean equations of motion 

< ...8S/8x... >= (12) 

The measure is given by 

»-\y) = (v^) D lim j exp(-%,*]/V^) (13) 



This expression is motivated by requiring that (i) the transfer matrix is the identity 
operator in the e — > limit; and (ii) symmetries of the action become symmetries of 
the integration measure. Finally, the evolution operator H in (||) is obtained by the 
e — ► limit 

H 



lim if f 



(14) 



followed by a rotation back to Lorentzian signature (5 — > rf) . 

This defines the quantization procedure. Next I will compute the operator H and 
obtain the classical limit of @, first with a flat space metric = rj^, then with a 
curved space metric. 

Let z^ = x M — y^. Since a rotation to Euclidean signature has been performed, we 
can easily evaluate the integral of equation (|i"3|) 



S[y,x] = 

H~\y) = J d D zexp[ 



m ^5 fll/ z' x z u = m \z\ 

mo\z\, 2n D / 2 



7h J T(D/2) 



{D-m—) D 

m 



and eq. 



(15) 



ip(y,T + e) 



J d D z fj, exp[ 
1 d 2 ip 



m \z\ 



+ 



eh { 

z fl z u + 0(z 3 



2 dy^dy 



dip 
dyv 



(16) 
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Identifying H t , and taking the limit and signature rotation prescribed in (0), we 
obtain 



H = -h 2 ^^r ] ^d,d u (17) 



and it is easy to see that the eigenstates of H 
are all solutions of the Klein-Gordon equation 



(FrTdA + = (19) 



where the mass-shell parameter 



D + r 



m 2 = (20) 
D + l K ' 



is proportional to the eigenvalue of the evolution operator H. 

The important point to notice, in this very simple example, is that the transfer 
matrix formalism has not destroyed the constraint, since each stationary state satisfies 
the usual Klein-Gordon equation. What is unusual is that the dimensionful parameter 
in the constraint, in this case m 2 , is different for each stationary state. In effect, m? 
has become a dynamical quantity, like energy, rather than a fixed parameter. 

The classical equations of motion can be obtained by WKB methods, or more 
simply by just making the replacement 



H[x»,-iKdr] -> H[af iPtl ] 

H c [x,p] = \im H[x,p] (21) 
h— >o 

and applying the classical equations 

H =£ dx^ = dH 1 dp }L = _dHc ^ 

dr dp^ ' dr dx^ 

where £ is fixed from the initial conditions. This gives us 

dr^ D + 1 

VTP0» = -M 2 , -T- = —±-tTp„ (23) 

dr mi 
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where M 2 = —1m\Ej (D + l). This is equivalent, of course, to the usual classical equa- 
tion of motion for a free relativistic particle moving in flat space, which is obtained 
by the rescaling 



ITLqS 



■ ~ M(D + 1) (24) 
where s is the proper time of the particle worldline. The proportionality of r to 
the proper time s has no special significance since, as noted below eq. the 
parametrization can be modified simply by choosing a different e(r) ^ 1. 

The non-relativistic limit of this theory is equally trivial, but still instructive. 
Making factors of c explicit, the general solution to the evolution equation fl5|) (with 
p° > 0) is 



J dm J dp f(m,p) exp |i[ ^ ~ 2 m 2 r — -^/m 2 c 4 +p 2 c 2 t + p ■ (25) 

where t = x°/c. For simplicity, suppose that: i) f(m,p) factorizes into f(m,p) = 
g(m)h(p); ii) g(m) is very sharply peaked around some value M; and iii) h(p) is 
negligible unless \p\ « Mc. In that case, again rescaling 



2^2 



M{D + 1) 

so that s has units of time, the solution can be written 



(26) 



ip(x^,s)^(p(s,t)if> NR (t,x) (27) 



where 



,m 2 c 2 



f(s,t) = I dm g(m)exp[i(^rs -mc t)/h] 



f P 
ip NR (t,x) = J dph{p)exp[-i{—t-p-x)/h] (28) 

and, of course, ipNR is a solution of the non-relativistic Schrodinger equation for a 
free particle of mass M. 

Now if g(m) and f(p) are both, e.g., gaussians, then the spreading of the wavepacket 
in both t and x is given by 

Ax(t) = [Ax 2 + (^) 2 ] 1/2 (29) 
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where 

AmAto ; , ApAxo ~ h (30) 

<r 

Since t =< t >= s, we can write A£(s) = At(t). Finally, it is consistent with the 
AmAto uncertainty relation to take both Am and Ato proportional to 1/c. This 
means that the dispersion of the wavepacket in the time direction 



Am = [Ati+(^m i/2 

~ -^0 (31) 

c 

goes to zero in the non-relativistic, c — > oo limit. 

The lesson of this simple exercise is that, in the c — > oo limit, the one- dimensional 
universe has acquired a clock. By a " clock" , I mean a non-stationary observable (or 
set of observables) whose dispersion is negligible, and whose evolution is independent 
of the other observables in the system. Measuring such an observable gives a value 
for the evolution parameter, which itself is not an observable. In the one-dimensional 
universe, the "clock" is the t-component of the field a; M (r); this is because its disper- 
sion can be made to vanish, and its value becomes perfectly correllated with r, in 
the non-relativistic limit. Equivalently, we can just say that the observable t behaves 
classically in the c — > oo limit. 

The extension of the transfer matrix formalism to a 1-dimensional universe with 
a curved field-space metric is straightforward. The reparametrization invariant 
action in this case is 




S = -m ldr X l-g^^- (32) 

Again we rotate from Lorentzian to Euclidean signature, extract an i, and obtain 
the action for Riemannian metrics. To evaluate the integrals in eq. (0) and (p~3|) , 
it is useful to introduce Riemann normal coordinates £ M around the point y^. The 
classical action for a trajectory running from £f = to £2 = £ is 



A^0,£]=m 0V /<Ue£ b + O(a (33) 

The 0(£ 5 ) terms will not contribute to H in the e — > limit, and can be dropped. 
The measure is then 
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det[^]«p[-mo|ei/V^] 



^■{D-l)^)**-* (34) 



F(D/2) ■ m ^y) ^y) 



The transfer matrix is computed as in the flat-space case, again with the help of 
Riemann normal coordinates 



ip(y,T + e) 



d A z 



a 



\lg{y + z) exp[-m ^g^z^z u / VeK]ifj(y + z, r) 
J (J o 

1 d 2 ip 



+ 



2 



(35) 



where R is the curvature scalar. Transforming back from Riemann normal coordinates 
and taking the limit fll4]) gives 



D 



" = -^ fi2 >^ + ^« (36) 

where g = |det(^ /Ji/ )|. The appearance of the curvature scalar term in the evolution 
operator H is related to the choice of measure, and, in turn, the ordering of operators 
in eq. (|36|) . There are other, more complicated choices of measure that could have 
been made in this problem, but these would only affect the coefficient of R (c.f. ref. 
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The classical equations of motion are again obtained either from the WKB ap- 
proximation, or else by the prescription (21), which gives 



~2mf 9 W 



(37) 



'o 

Setting H c = £, and M 2 = —2m^£/{D + 1), and also rescaling r according to 
we recover the the classical equations of a relativistic particle moving in curved space, 
with the identification of the evolution parameter r as the proper time. 

At the quantum level, eigenstates of the mass-shell parameter are stationary states 
of the evolution operator H . The remarks made above for the flat-space case, noting 
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that all stationary states obey the constraint equation but with different values of 
the mass-shell parameter, of course apply to the curved space example also. Non- 
stationarity, in the one-dimensional universe, must be attributed to dispersion in the 
mass-shell parameter M 2 . 

As a final remark we may ask whether, since the evolution parameter is identified 
with proper-time in the classical limit (given the choice e(r) = 1 in eq. (]?])), the 
proper-time could have been used as an evolution parameter from the beginning. This 
might be done by replacing S[y, x] and e in eq. @ with a weight S[(y, s + As), (x, s)], 
representing the minimal action between points x and y of paths constrained to 
have proper time As. The problem with this is that the amplitude between (x,si) 
and (y, S2) would not, in general, be dominated by the classical path, unless the 
proper time difference S2 — S\ happened to correspond to the path length of the 
geodesic between x and y. The problem can be fixed, for Green's functions, by 
integrating over S2 — si (see, e.g., @]), but then the proper-time loses its function as 
an evolution parameter. In higher dimensions the situation is further complicated by 
the fact that a proper-time slicing of a simply connected Riemannian D- manifold into 
(D-l)-manifolds can introduce spurious singularities. Under proper-time evolution, 
a simply connected (D-l)-manifold will in general evolve into a set of disconnected 
(D-l)-manifolds ||. It seems unlikely that such an approach would yield a hermitian 
evolution operator, although it may have other applications (c.f. |§). 



3 Minisuperspace Models 

Next we consider actions of the form 



S = Jdt( Pn ^£-NH) 



U = -^-G nm p nPm + m V(q) (38) 

where m is some dimensionful parameter, and the supermetric G nm has Lorentzian 
( — h+...+) signature. Actions of this kind typically arise in minisuperspace models of 
quantum cosmology. To compute S[q + Aq, q], begin with the Hamiltonian equation 

q n = N ^L = ^-Q nm Vm (39) 
opn m 
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insert this into the constraint equation 

U = ^G nm q n q m + m V = (40) 

and solve for the lapse 

N=[-^rjG nm q n q m } 1/2 (41) 



2V 



Then 



S[q,q + A] = [ dt^G nm q n q m 
Jo JM 

= -m J dt^J-2VG nm q n q r 



-m J-2VG nm Aq n Aq T 



= -m ^/-g nm Aq n Aq m (42) 
where we define a modified supermetric 

Qnm = 2VG nm (43) 

From here, the procedure follows exactly the same steps as in the curved metric 
example of the previous section. Again, rotate the signature of the supermetric G nm 
(not the spacetime metric g^) to Euclidean signature and extract a factor of % to 
obtain the "Riemannian" action. Then, introducing Riemann normal coordinates 
in minisuperspace around the point q + Aq, evaluate the relevant integrals for the 
measure and transfer matrix. The final result for the evolution operator is 

2ml VGdq n dq m 6m 2 K ' 

where, in this case, D is the dimensionality of minisuperspace, 1Z is the curvature 
scalar formed from the modified supermetric ([i"3"D, and 



^ = |det(^ mn )| (45) 

It should be noted that H is hermitian in the measure /z(g) = \[Q. This is de- 
spite the fact that, in quantum cosmology models, V(q) is not positive definite, and 
therefore S[q,q + Aq] can be imaginary, even after rotation of G nm . Nevertheless, 
carrying out the relevant integrals and rotating back to Lorentzian signature, one 
still finds that H is hermitian in the appropriate measure. This can be understood 
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from the fact that, for modified supermetrics Q mn of Euclidean signature, the operator 
H e is hermitian by construction, since S[q, q'\ is real and symmetric. But since the 
hermiticity of H e does not depend on the precise functional form of Q m n, the contin- 
uation to arbitrary signature can only upset hermiticity if it introduces factors of i in 
H . These factors of i are avoided by making the standard Euclidean — > Lorentzian 



continuation Jdet(g) — > J-det(g), in the one-dimensional Universe example, and 



cfe£((? mri ,) — > y\det(Gmn)\, m the minisuperspace case, for both H and fx. 
The Schrodinger equation for stationary states is now 



^ T ^w^ e '' m w + ^ n ^ £ * (46) 



or, equivalent ly 



1 V 8 i 1 8 1 1m 2 F 

-h 2 ^=^J\VG\^-G nm ^- + -h 2 VTZ - ^-V]* = (47) 

2 x f\VG\dq nyi { V dq rn 3 D + 1 K ' 



This equation may look more recognizable if we identify 

M 2 = (48) 
D+ 1 v ; 

Then reorder the factors in the Laplacian and drop the curvature scalar term to obtain 

" ^t^^W + MV] * = ° " (49) 
where the quotation marks mean that this equation is correct up to operator-ordering 
terms. From this it is easy to see that the equation satisfied by stationary states, 
eq. (0), is simply the Wheeler-DeWitt constraint equation, with the dimensionful 
parameter identified as M 2 = — 2m^£ /{D+l), and with a particular choice of operator 
ordering. Once again, we see that all stationary states satisfy a Wheeler-DeWitt 
equation, but with different mass parameters. 

As in the last section, the classical equations of motion are obtained from the 
classical Hamiltonian given by (^l|): 



He = T-^G nm PnPm (50) 



D + l 

o ' 

and, setting H c = £ and M 2 = —2rn^£/{D + 1), we obtain the classical constraint 
equation 

H M = ^G nm PnPm -MV = (51) 
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Hamilton's equations give us 

dH r D 



dq n 
dr 

dp n 
dr 



1 



dp n 2m\V 
dH r D 



Qnm 
+ 1 



dq n 

(D + 1)M 
ImlV 



Ami 



1 dGv 
V~dq 

dG ij 



PiPj 



2M dq 



M 



1 dV 
V^dq^ 
dV_ 

dq n 



G ij piPj 



(52) 



which can be rewritten as 



H M = 

dq n , s &Hm 

dr dp n 

~d7 = (53) 

where 

. (D + 1)M , . 

m = (54) 

and where H-m, given in (|5l|) , is the minisuperspace Hamiltonian 7i with the param- 
eter mo replaced by M. 

Apart from the factor of f(q), equations ( |53| ) are simply the classical equations 
of motion of the minisuperspace action (with m — > M). So long as V is non-zero, 
f(q), like the lapse function N(t), is irrelevant in determining the classical trajectory, 
which depends only on the directions of the vectors d T q n and d T p n in phase space, and 
not on their magnitudes. The magnitudes of these vectors only determine the rate 
(compared to some analog of proper time in minisuperspace) at which the evolution 
parameter runs along the classical trajectory. 

This point can be made in another way. Suppose we pick an initial point in phase 
space {q,p}o and solve the equations (|5^). Denote the solutions g"(r) and p n (r). 
Then, define 

N(r) {q , p}0 = f(q(r)) (55) 

It is easy to see that that ^"(r) and p n {'T) are a solution of Hamilton's equations for 
the original minisuperspace action ([38]) 

Hm = 
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dq n &H M 

d P n ATI \ 31~~Lm fr . a \ 



with mo replaced by M, and with a particular choice (|55| ) for the lapse function 
(which will be different for each classical trajectory). 

4 Quantum Gravity 

The action of quantum gravity, in the ADM decomposition, is 



S = <Tx 







ik. 



G^ki — 2^j=(9ik9ji + 9iWjk — 9ij9ki) (57) 

where ^fg is the root determinant of the three- metric gij. At the classical level, 
k,q = 167tGn, where Gat is Newton's constant. It will be seen that this identification 
is modified at the quantum level, much as wig was replaced by M 2 in the examples of 
the preceding sections. 

The presence of the shift functions iVj is a serious complication, as compared to 
the minisuperspace models of the preceding section, where the shift functions were 
absent. The problem is that the classical equation of motion 



dgij 



2K 2 NG ljnmP nm + N i;j + N j;i (5* 



dt 

which are used to solve for the momenta in terms of dtgij, contain derivatives of the 
iVj. These shift functions can be solved for, in principle, by substituting the expression 
for obtained from (|5^) into the supermomentum constraints 7i l x = 0, which gives 
the iVj in terms of the lapse N, and then substituting these iVj into the Hamiltonian 
constraint 7i x = to solve for the lapse. However, since the iVj are determined, in this 
way, by complicated partial differential equations, this procedure leads to intractable 
expressions, and there is no simple form for S\g', g]. 
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A great simplification is achieved if we instead set iVj = from the beginning. In 
that case the supermomentum constraint is not obtained by extremizing the action, 
and must somehow be recovered by imposing an operator constraint on physical states 

Qx\Pij,gn]v = o (59) 

where the subscript x indicates that there is a separate constraint at each point, and of 
course these must be mutually consistent, as well as consistent with the r-evolution 
operator. These constraints will be obtained below. For the moment, we just set 
Ni = and proceed. Denote 



{a = 1 - 6} <-> i < j} 

q a (x) <-» gij (x) 

Pa[X) \2pv{x) (i<j) 

G ab (x) <- G ijnm (x) 



(60) 



in order to compress the number of indices somewhat. Hamilton's equations with 
Ni — give us 



1 



Pa 



2k 2 N 



Inserting this into the constraint equation 

1 



= ^ 



G ab q b 



1 



x 



where 



U = -R + 2A 



G ab q a q b 



1/2 



and solving for the lapse, gives 

N = 
so we have 

AS = S[q',q} = -\ I d 3 x [ At dt^-g- 2 UG ab qOq b 
Kq J Jo 



= -\ I d 3 x^-(g-2UG ab ) Aq*Aqb + 0(Aq 2 ) 
k J 

_J_ j d 3 x(^) (^-(g ab ) Aq"Aq» + O(Aq 2 ) 



(61) 

(62) 
(63) 

(64) 



(65) 
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where Aq a = q a — q' a , and we define 

g ab = ^-UG ah (66) 

The notation (..) means that the quantity in parenthesis is to be evaluated at Aq = 0; 
i.e. at q'. 

Once again, we implicitly rotate the signature of the superspace metric G a b to 
Euclidean signature, and extract a factor of i, before evaluating the transfer matrix. 
Note that it is the signature of superspace, not the signature of spacetime, which 
is rotated to the Euclidean value. At the end, of course, we rotate back to the 
usual signature of superspace, according to the prescription ([L4]). It was found, in 
the minisuperspace example of the previous section, that the evolution operator is 
Hermitian in an appropriate measure fi(q), despite the fact that the potential term 
V was not positive definite. It is expected, for similar reasons, that hermiticity of H 
will also be maintained in full quantum gravity, where U is not positive definite, f] 

The transfer matrix is obtained from 



i)(q',T + e) 



+ 



Dqfi(q)e- AS/VIK 
- J d 3 xd 3 y 



ip(q') + I d 3 x 



5ip 



Aq a {x) 



2J~ ~ a \8q a (x)8q b (y) / 
l + etTo + H + Tsl + Ote 2 ) 



Sq a (x 
Aq a {x)Aq\y) + ... 



(67) 



where the T n represent terms with n derivatives of if). To find these terms, we need 
to evaluate 



< Aq a (x 1 )Aq b (x 2 ) > = J D(Aq) (fi) Aq a (x 1 )Aq b (x 2 ) 



x exp 



~ / d 3 x{^) ( ^(g ab ) Q Aq a Aq b / Vrfi 



(68) 



This quantity is highly singular, and, at this point, regularization is unavoidable. 

Unfortunately, the problem of regularizing a functional integral 
non-perturbatively, in such a way as to preserve diffeomorphism invariance, is un- 
solved. So I will have to resort to a naive regulator, and replace the continuum 



6 An actual proof of hermiticity, however, would involve the regularization/operator-ordering is- 
sues discussed below. 
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degrees of freedom, labeled by x, by a finite set, labeled {n}, which are taken to 
represent regions of equal volume. We make the correspondences 



Aq a (x) 
5 

Sq a (x) 
Dq 



Aq a (n) 
V Np 

— Y 

JV P n=l 



V dq a (n) 
Ed D q( 



n 



(69) 



where V is the three-volume. For pure quantum gravity, there are six independent 
components of <?y, so D=6. Obviously, with this naive discretization we lose diffeo- 
morphism invariance of the integration measure, and we can expect to make mistakes, 
even after taking the continuum limit, on certain expressions which depend crucially 
on the invariance of the measure. In the transfer matrix formulation, it is the operator 
ordering, and the presence of terms analogous to R in eq. (^) (or 1Z in ([|7])), which 
are sensitive to the measure. Operator-ordering problems are notorious in quantum 
gravity, and I will not try to solve them here. From here on, only the principal term 



of the evolution operator (coming from T2 in eq. (|67|) ), will be determined explicitly. 
This term is relatively insensitive to the measure, and is the only term which is im- 
portant for the classical limit. But even for the T 2 term, the ultimate justification for 
the prescription (|69"D above will be a posteriori. 
With the above caveats duly noted, we find 



< Aq a (n)Aq\m) > = J J[ d D q{k) (fi) Aq a (n) Aq b 



m) 



cxp 



V 



v 2 E \lG ab )zAq«{k)Aq\k)/^K 



K 4 N 2 



then 



d 2 ^ 



€T2 2^^~dq-(n)dq\m) 



(34 



K 4 N 2 



v 2 

< Aq a (n)Aq h (m) > 

d 2 ij 



(70) 



'06 



eh 



V 



V 2 V dq a {n)dq b {n) 



d 6 x U- L G 



— l/^ab 



dq a dq b 



(71) 
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where (3 is a divergent, dimensionless constant 

fi = l(D + l)N p (72) 
proportional to the number of degrees of freedom. Finally, from Ti we get 

H = -h 2 ^ [ d 3 x U-^-^— (73) 
V J 5q a 5q b v ' 

where, again, operator-ordering/measure terms have been dropped. 

However, this evolution operator does not yet give us general relativity, not even 

in the classical limit. This is easy to see from, e.g., the WKB approach, since the 

Hamilton- Jacobi equation corresponding to ([73]) is 



I y J dXU Cr s ^ 5 ^ {(4) 

which is not the Einstein-Hamilton- Jacobi equation of gravitation. Now, although the 



transfer matrix formulation as presented in eq. @ and (|T^) was guaranteed to recover 
the correct classical limit, it is not hard to see what went wrong in this case. The 
difficulty, as already mentioned above, comes from setting the shift functions N{ = 
at the beginning; this means that the supermomentum constraints TC X = were lost 
from the start. One could try to simply reintroduce these as physical state constraints 
TC l x ^ = 0, but this still doesn't lead to the Einstein-Hamilton- Jacobi equation. 

Note however, that because of the factor 1/V in the r-evolution operator H of 



eq. (|73f) , the Schrodinger equation can be written in the form 

d 3 xQ x ^ = (75) 



where 

Q x = -h^p-^^ - ih^ (76) 

The extra constraints which need to be imposed on the physical states, which then 
generate the usual constraint algebra of general relativity, are simply 

Q X V = (77) 

at every point x. It is not hard to see why. Consider an arbitrary solution of the 
Schrodinger equation 

*(q,T) = J2*e e ~ ieT/h M<l) (78) 
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Since the a £ are arbitrary, the Q-constraint (|77|) requires that for each stationary 
state 

l- h2Gab J^7 b + V~9 u)<l> £ =0 (79) 



Sq a 5q b v 



Identifying 



K 2 = 16ttG 



N 




p = !/-/>-*- (80) 

as Newton's constant and the Planck mass, respectively, we see that eq. ( |79|) is just 
the Wheeler-DeWitt equation 

H x (^e = {-^ 2 G ab ^ b + ±V9U] *e = (81) 

with Newton's constant inversely proportional to \/—E. 

Finally, we invoke the Moncrief-Teitelboim interconnection theorem ||10| , which 
says that if a state satisfies the (Wheeler-DeWitt) Hamiltonian constraint fl8lD at 
every point x, then that state also satisfies the supermomentum constaints 

W x <&s = (82) 

at every point x. In this way, the supermomentum constraints that were lost at the 
outset by fixing iVj = have been recovered. Further, given that the Hamiltonian 
and supermomentum constraints are consistent (commutators close on the Poincare 
algebra), and that the Q x constraints ([77|) are implied by the Hamiltonian constraints 
(pl|), it follows that the Q x constraints are consistent not only with the evolution 
operator H, but also with each other. 

The Einstein-Hamilton- Jacobi equations follow directly from a WKB treatment 
of the Wheeler-DeWitt constraints, and the classical limit, obtained by replacing 
—ih5/5q a with p a , follows in complete analogy to the minisuperspace case. Since the 
proper constraint algebra has been obtained, it is fairly obvious that that the correct 
classical equations must fall out. But it is still nice to see this explicitly. We begin 
with the classical limit of the evolution operator 



H c = - J d 3 x ^4U^G ab PaPb 
= — I d 3 x II 



J3) 
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Begining from a set of initial data {q,p}o = {9ij,P^}o consistent with the supermo 
mentum constraints 

n i [{q, P u = o 



we have 

S = H c [{q,p} } 

and consistency of the initial data with the Q x = constraints gives us also 

Hex = \/g£ 

Then Hamilton's equations, derived from H c , are 

dq a (x) (34 



(84) 
(85) 

(86) 



2G ab p b 



(It 



uv 



3 (3 4 5 2 , 
dX UV^5p a (x) Hx {K ' 



and 



dPa(x) 

dr 



1 5V 



+ - I d 3 x 



V 2 5q a (x 



d 3 x'H cx , / d 3 x 



, 1 SU oM 



V J U 2 5q a (x 



P4G bc PbPc 



G bc PbPc 



Applying eq. fl8"6|) 

dp a {x) 
dr 



U 5q a (x] 



- / d 6 x 



UV 1 5q a (x 



+(3k1 



S 



G bc PbPc 



x 



d 3 x 



, P4 S 



UVk 2 Sq a (x 



Sq a (x) 



K 2 G bc p b p c + -^U 



(89) 



The evolution parameter r has units of mass x length. To give it conventional units, 
rescale r — ► t = r/Mp, and define 



F x [q] = 



M p [3k 4 
U{x)Yk< 



(90) 
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Then the full set of gravitational field equations derived from H c plus the con- 
straints is 



H x (k 2 ) = 

K = (91) 

To see that this is classical Einstein gravity, let q a (x l , t) and p a {x l , t) be a solution of 
(pl|) for some set of initial data {q,p}o = {gij,P^}o compatible with the constraints. 
Then, as in the minisuperspace case, define 

N(x,t) {p , q}o =F x [q(t)} (92) 
and one sees that q a (x\t) and p a (x\t) is a solution of 

dq a {x) SH E instein 

dt Sp a (x) 

dp a (%) _ SHEinstein 

dt 5q a (x) 

H Emstem = Jd 3 xN(x,t) {M}o H x (K 2 ) 

H x (k 2 ) = 

K = (93) 



which are simply the Einstein field equations in Hamiltonian form, with lapse (]9 
which is in general different for each classical solution, and shift Ni = 0. 



5 Discussion 



To summarize: given an initial " wavefunction of the Universe" *&[q, To], it has been 
shown how to trace its subsequent evolution along a trajectory in the Hilbert space 
of physical states, with the trajectory parametrized by r. The r-evolution equation 
for quantum gravity, up to factor-ordering terms, is 

ihd T V[q,T] = H^[q,r] 

= ^Jd-xU-^^-^A (94) 
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where the space of physical states is spanned by the eigenstates 



H$ £ [q]=E$ e [q] 



(95) 



each of which satisfies the constraint algebra 



H x (k 2 )$ s 

K 4 












(96) 



and where H. x {k 2 ) denotes the superHamiltonian constraint operator, with Newton's 
constant Gn = k 2 /167t. Modulo operator-ordering issues, every stationary physical 
state in the transfer matrix formulation satisfies the usual constraint algebra of general 
relativity, but each with a different value of Newton's constant, where G 2 N is inversely 
proportional to the eigenvalue £ of the evolution operator. The number of degrees of 
freedom, in this formulation, is therefore 2 x oo 3 + 1; i.e. two degrees of freedom per 
point, which is the degrees of freedom of states satisfying the constraint algebra (|96|). 
and one extra degree of freedom corresponding to the Planck mass. This is only one 
more degree of freedom (overall, not per point) than in the standard formulations; 
there should be no danger of, e.g., the graviton aquiring a mass. 

Thus, in the transfer-matrix formulation, time-evolution of states is recovered at 
a modest price: the Planck mass (inverse Newton's constant) becomes a dynamical 
quantity, analogous to energy in non-relativistic quantum mechanics. As is the case 
for the time parameter in non-relativistic quantum mechanics, or in quantum field 
theory on a fixed spacetime background, the evolution parameter r is only an evolution 
parameter; it is not an observable, and there is no operator acting on the Hilbert space 
which corresponds to r. This avoids the problems encountered in relativistic quantum 
mechanics, as well in standard formulations of canonical quantum gravity, where the 
approach is to identify one of the operators in the theory as an evolution parameter. 

We note that the spectrum of the evolution operator H in fl9"3| ) is unbounded from 
below; this is true of all the systems considered in this paper. For closed systems 
(such as the Universe) which do not interact with anything external, this absence of 
a ground state is not a problem. The eigenvalue £ is a constant of motion, and its 
distribution cannot change. 

In non-relativistic quantum mechanics, any non-stationary state has a certain dis- 
persion in its energy. For quantum gravity, the corresponding statement is that since 
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the Universe is non-stationary, there must be a certain dispersion in S, the eigenvalue 
of the evolution operator. This implies dispersion in the Planck mass Mp, or, equiv- 
alently, Newton's constant. Since Mp oc — £ is conserved by the evolution equation, 
the fractional dispersion of Newton's constant AGn/Gn should be a dimensionless 
constant of nature. Depending on how large this dispersion is, there could conceivably 
be observational consequences. I hope to return to this question at a later time. 
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